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Abstract 

The linear theory of discrete Riemann surfaces is applied to polyhe- 
dral surfaces embedded in M^. As an application we compute the period 
matrices of some classical examples from the surface theory, in particular 
the Wente torus and the Lawson surface. 



1 Introduction 

Finding a conformal parameterization for a surface and computing its period 
matrix is a classical problem which is useful in a lot of contexts, from statistical 
mechanics to computer graphics. 

The 2D-Ising model 119j|Bil9] for example takes place on a cellular decomposi- 
tion of a surface whose edges are decorated by interaction constants, understood 
as a discrete conformal structure. In certain configurations, called critical tem- 
perature, the model exhibits conformal invariance properties in the thermody- 
namical limit and certain statistical expectations become discrete holomorphic 
at the finite level. The computation of the period matrix of higher genus sur- 
faces built from the rectangular and triangular lattices from discrete Riemann 
theory has been addressed in the cited papers by Costa-Santos and McCoy. 

Global conformal parameterization of a surface is important in computer 
graphics [TBI HH HI HHl EZl HZ] in issues such as texture mapping of a flat picture 
onto a curved surface in R'^ . When the texture is recognized by the user as a nat- 
ural texture known as featuring round grains, these features should be preserved 
when mapped on the surface, mainly because any shear of circles into ellipses is 
going to be wrongly interpreted as suggesting depth increase. Characterizing a 
surface by a few numbers is as well a desired feature in computer graphics, for 
problems like pattern recognition. Computing numerically the period matrix of 
a surface has been addressed in the cited papers by Gu and Yau. 
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This paper uses the general framework of discrete Riemann surfaces the- 
ory [Ml [131 mi d] and the computation of period matrices within this frame- 
work (based on theorems and not only numerical analogies). Wc translate these 
theorems into algorithms. Some tests arc performed to check the validity of the 
approach. 

We start with some surfaces with known period matrices and compute nu- 
merically their discrete period matrices, at different level of refinement. In 
particular, some genus two surfaces made out of squares and the Wente torus 
are considered. We observe numerically good convergence properties. Moreover, 
we compute the yet unknown period matrix of the Lawson surface, recognized 
it numerically as one of the tested surfaces, which allowed us to conjecture their 
conformal equivalence, and finally prove it. 

2 Discrete conformal structure 

Consider a polyhedral surface in M.^. It has a unique Dclaunay tesselation, 
generically a triangulation [5]. That is to say each face is associated with a 
circumcircle drawn on the surface and this disk contains no other vertices than 
the ones on its boundary. Let's call F the graph of this cellular decomposition. 
To its vertices, Fi its edges and complete it into a cellular decomposition with 
F2 the set of triangles. Each edge {x,x') = e g Fi is adjacent to a pair of 
triangles, associated with two circumcenters y,y'. The ratio of the (intrinsic) 
distances between the circumcenters and the length of the (orthogonal) edge e 
is called p{e). 



Following [19] , we call this data of a graph F, whose edges are equipped with a 
positive real number a discrete conformal structure. A discrete Riemann surface 
is a conformal equivalence class of surfaces with the same discrete conformal 
structure. It leads to a theory of discrete Riemann surfaces and discrete analytic 
functions, developed in [THinilinKinillllllllll], that shares a lot of features with 
the continuous theory and these features are recovered in a proper refinement 
limit. We are going to summarize these results. 

In our examples, the extrinsic triangulations are Delaunay. That is to say the 
triangulations come from the embedding in K."^ and the edges {x, x') G Fi of the 
triangulation are the edges of the polyhedral surface in M'^. On the contrary, 
the geodesic connecting the circumcenters y and y' on the surface is not an 
interval of a straight line and its length is generically greater than the distance 
IIj/ — y'll in R^. The later gives a more naive definition of p, easier to compute 
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that we will call extrinsinc, because contrarily to the intrinsic p, they depend 
on the embedding of the surface in MP and are not preserved by isometrics. 
For surfaces which are refined and flat enough, the difference between extrinsic 
and intrinsic distances is not large. Wc compared numerically the two ways 
to compute p. The conclusion is that, in the examples we tested, the intrinsic 
distance is marginally better, see Sec. 14.21 

The circumcenters and their adjacencies define a 3-valent abstract (locally 
planar) graph, dual to the graph of the surface, that we call F*, with vertices 
Fq = F2, edges F* ~ Fi. We equip the edge {y,y') = e* G F*, dual to the 
primal edge e G Fi, with the positive real constant p{e*) = l/p{e). We define 
A := r©F* the double graph, with vertices Aq = FqUFq and edges Ai = FiUF*. 
Each pair of dual edges e, e* e Ai, e = {x, x') G Fi, e* = (y, y') e F*, are seen 
as the diagonals of a quadrilateral {x, y, x', y'), composing a quad-graph 0, with 
vertices Oo = Aq, edges Oi composed of couples {x,y) and faces O2 composed 
of quadrilaterals (x, y, x' , y'). 

The Hodge star, which in the continuous theory is defined by *(/ dx+g dy) = 
~g dx + f dy, is in the discrete case the duality transformation multiplied by 
the conformal structure: 

*a := p{e) j a (1) 



e* 



A 1-form a E C^{A) is of type (1,0) if and only if, for each quadrilateral 



{x,y,x',y') G O2, I(yy)(^ = ^Pi^'X') I[x,x')^' *° if 



*a = 



Similarly forms of type (0, 1) are defined by *a — +ia. A form is holomorphic, 
resp. anti-holomorphic, if it is closed and of type (1,0), resp. of type (0, 1). A 
function / : Aq ^ C is holomorphic iff d\f is. 

We define a wedge product for 1-forms living whether on edges Oi or on 
their diagonals Ai, as a 2-form living on faces O2. The formula for the latter is: 

«A/3:=i I I a I P- I a I (3\ (2) 



2 

(a:,J/,x',i/') \lx,x') (y,y') (y,y') (x.x') 

The exterior derivative d is a derivation for the wedge product, for functions 
/, g and a 1-form a: 

d{fg) = fdg + gdf, d{fa) ^ dj A a + f da. 

Together with the Hodge star, they give rise, in the compact case, to the usual 
scalar product on 1-forms: 

{a,(i):^([ a A *^ = (*a, */3) = (A^ = i ^ p(e) / a / /3 (3) 



eeAi 



The adjoint d* = — * d * of the coboundary d allows to define the discrete 
Laplacian A = d* d -I- d d* , whose kernel are the harmonic forms and fimctions. 
It reads, for a function at a vertex a: G Aq with neighbours x' ^ x: 



(A/)(x)= ^p(x,x')(/(x)-/(x')). 
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Hodge theorem: The two ±z-eigenspaces decompose the space of 1-forms, 
especially the space of harmonic forms, into an orthogonal direct sum. Types 
are interchanged by conjugation: a S C^^'^^{A) a G C^'^'^\A) therefore 

{a,f3) = (7r(i.o)a,7r(i^o)/3) + (7r(oa)a, 7i'(o,i)/5) 
where the projections on (1, 0) and (0, 1) spaces arc 

7i"(i,o) = ^(Id + «*), 7r(o,i) = ^{Id-i*). 

The harmonic forms of type (1,0) are the holomorphic forms, the harmonic 
forms of type (0, 1) are the anti-holomorphic forms. 

The norm of the 1-form df, called the Dirichlet energy of the function /, 
is the average of the usual Dirichlet energies on each independent graph 

En{f):=\\dfr^idf,df) = ^ J2 Pix,x')\f{x')-f{x)f (4) 

{x,x'}eAi 
_ Epjflr) + Eoiflr^) 
2 

The conformal energy of a map measures its conformality defect, relating these 
two harmonic functions. A conformal map fulfills the Cauchy-Riemann equation 

*df = -idf. (5) 

Therefore a quadratic energy whose null functions are the holomorphic ones is 

Ec{f):^^\\df^i*dff. (6) 

It is related to the Dirichlet energy through the same formula as in the contin- 
uous case: 

Ecif) ^l{df-i*df,df-i* df) 

= ^\\dff + ^J\-i*dff+Rc{df, -i^dj) 

= ||d/||2+ Im // d/Arff 

J J()2 

= Eoif)^2A{f) (7) 



where the area of the image of the application / in the complex plane has the 
same formulae (the second one meaningful on a simply connected domain) 

^if)='oll dfhd!=-i fW-Jdf (8) 

as in the continuous case. For a face (x, y, x' , y') G O2, the algebraic area of the 
oriented quadrilateral ( /(x), /(x'), /(y), f{y') ) is given by 



ix,y,x',y') 



dfAdf^z Im ((/(x') - /(.T))(/(y') - fiy))) 
= -2iA(f{x)J{x'),f{y)J{y')] 
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When a holomorphic reference map z : Aq ^ C is chosen, a holomorphic 
(resp. anti-holomorphic) 1-form df is, locally on each pair of dual diagonals, 
proportional to dz, resp. dz, so that the decomposition of the exterior derivative 
into holomorphic and anti-holomorphic parts yields dfAdf = + dzA 

dz where the derivatives naturally live on faces and are not be confused with 
the boundary operator d. 

3 Algorithm 

The theory described above is straightforward to implement. The most sensitive 
part is based on a minimizcr procedure which finds the minimum of the Dirichlet 
energy for a discrete Riemann surface, given some boundary conditions. Here is 
the crude algorithm that we are going to detail. A normalized homotopy basis 
H of a quad-graph 0{S) is a set of loops G kcr9 

Basis of holomorphic forms (a discrete Riemann surface S) 

find a normalized homotopy basis H of (}{S) 

for all Hfc do 
compute and 

compute the real discrete harmonic form ujk on T s.t. = 7 o 

(check LOk is harmonic on T) 
compute the form *cjfc on T* 
(check ^LOk is harmonic on T*) 
compute its holonomies (/^r* *^k)k,i on T* 
do likewise for lo^ on F* 
end for 

do some linear algebra (R is a rectangular complex matrix) to get the basis 
of holomorphic forms {Ck)k = ^(Id + i *){uJk)k s.t. (/^r Cfc) = 4,^ 
define the period matrix Ilk,i (/^r* Cfc) 
do hkewise for (C^)fc and n^_^ := (/^r Q) 

Finding a normalized homotopy basis of a connected cellular decomposition 
is performed by several well known algorithms. The way we did it is to select a 
root vertex and grow from there a spanning tree, by computing the vertices at 
combinatorial distance d from the root and linking each one of them to a unique 
vertex at distance d — 1, already in the tree. Repeat until no vertices are left. 

Then we inflate this tree into a polygonal fundamental domain by adding 
faces one by one to the domain, keeping it simply connected: We recursively 
add all the faces which have only one edge not in the domain. We stop when 
all the remaining faces have at least two edges not in the domain. 



5 



Period Matrices of Polyhedral Surfaces 



Bobenko, Mercat & Schmies 




Then we pick one edge (one of the closest to the root) as defining the first 
element of our homotopy basis: adding this edge to the fundamental domain 
yields a non simply connected cellular decomposition and the spanning tree 
gives us a rooted cycle of this homotopy type going down the tree to the root. 
It is (one of) the combinatorially shortest in its (rooted) homology class. We 
add faces recursively in a similar way until we can no go further, we then choose 
a new homotopy basis element, and so on until every face is closed. At the end 
we have a homotopy basis. We compute later on the intersection numbers in 
order to normalize it. 




We compute the unique real harmonic form rj associated with each cycle H 
such that 77 = 7 o H. This is done by a minimizing procedure which finds the 
unique harmonic function / on the graph F, split along H, whose vertices are 
duplicated, which is zero at the root and increases by one when going across H. 
This is done by linking the values at the duplicated vertices, in effect yielding a 
harmonic function on the universal cover of F. The harmonic 1-form d f doesn't 
depend on the chosen root nor on the representative H in its homology class. 

4 Numerics 

We begun with testing discrete surfaces of known moduli in order to investigate 
the quality of the numerics and the robustness of the method. We purposely 
chose to stick with raw double 15-digits numbers and a linear algebra library 
which is fast but not particularly accurate. In order to be able to compare 
period matrices, we used a Siegel reduction algorithm |10| to map them by a 
modular transformation to the same canonical form. 



6 



Period Matrices of Polyhedral Surfaces 



Bobenko, Mercat & Schmies 



4.1 Surfaces tiled by squares 

Robert Silhol supplied us with sets of surfaces tiled by squares for which the 
period matrices are known [3 El [23] . There are translation and half- 
translation surfaces: In these surfaces, each horizontal side is glued to a hori- 
zontal side, a vertical to a vertical, and the identification between edges of the 
fundamental polygon are translations for translation surfaces and translations 
followed by a half-twist for half-translations. The discrete conformal structure 
for these surfaces is very simple: the combinatorics is given by the gluing con- 
ditions and the conformal parameter p = 1 is constant. 

The genus one examples are not interesting because this 1-form is then the 
unique holomorphic form and there is nothing to compute (the algorithm does 
give back this known result). Genus 2 examples are non trivial because a second 
holomorphic form has to be computed. 

The translation surfaces are particularly adapted because the discrete 1-form 
read off the picture is already a discrete holomorphic form. Therefore the com- 
putations are accurate even for a small number of squares. Finer squares only 
blur the result with numerical noise. For half-translation surfaces it is not the 
case, a continuous limit has to be taken in order to get a better approximation. 



Surface & Period Matrix 



Numerical Analysis 
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^/^vcrticcs 
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14 
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62 
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4096 
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^^vcrticcs 



10, 



22 



3.40 



O3 = ^ 



2 -1 
73 I -1 2 



94 



9.51 



382 



2.44 



1534 



6.12 



6142 



1.53 



Using 15 digits numbers, the theoretical numerical accuracy is limited to 8 
digits because our energy is quadratic therefore half of the digits are lost. Using 
an arbitrary precision toolbox or Cholesky decomposition in order to solve the 
linear system would allow for better results but it is not the point here. 
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4.2 Wente torus 




Grid : 40 X 40 Grid : 80 x 80 

Figure 1: Regular Delaunay triangulations of the Wente torus 

For a first test of the numerics on a an immersed surface in IR'^ our choice is 
the famous CMC-torus discovered by Wente [28] for which an explicit immersion 
formula exists in terms of theta functions [3] . The modulus of the rhombic Wente 
torus can be read from the immersion formula: 

T„ « 0.41300 . . . + 0.91073 ...ifn exp(a. 145045 ....). 

We compute several regular discretization of the Wente torus (Fig. [T]) and 
generate discrete conformal structures using pcx that are imposed by the ex- 
trinsic Euclidean metric of R"^ as well as pin which are given by the intrinsic 
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flat metric of the surface. For a sequence of finer discretizations of a smooth 
immersion, the two sets of numbers come closer and closer. For these discrete 
conformal structures we compute again the moduli which we denote by Tex and 
Tim and compare them with t^ from above: 



Grid 


l|Ti„ - T^ll 


i '^cx '^w 1 


10 X 10 


5.69- 10^^ 


5.00-10-^ 


20 X 20 


2.00-10-^ 


5.93-10-^ 


40 X 40 


5.11 • 10--* 


1.85 • 10-^ 


80 X 80 


2.41 • 10-* 


6.00-10-'' 



For the lowest resolution the accuracy of t^x is slightly better then the one 
of Tin. For all other the discrete conformal structures with the intrinsically 
generated pin yields significant higher accuracy. 

4.3 Lawson surface 





1162 vertices 



2498 vertices 



Figure 2: Delaunay triangulations of the Lawson surface 



Finally we apply our method to compute the period matrix of Lawson's 
genus 2 minimal surface in §^ [THj. Konrad Polthicr [52] supplied us with sev- 
eral resolution of the surface which are generated by a coarsening and mesh 
beautifying process of a very fine approximation of the Lawson surface (Fig. [2]). 
Our numerical analysis gives evidence that the period matrix of the Lawson 
surface is 

which equals the period matrix fl^ of the third example from Sec. 14.11 Once 
conjectured that these two surfaces are conformally equivalent, it is a matter 
of checking that the symmetry group of the Lawson genus two surface yields 
indeed this period matrix, which was done, without prior connection, in [1]. An 
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explicit conformal mapping of the surfaces can be found manually: The genus 2 
Lawson surface exhibits by construction four points with an order six symmetry 
and six points of order four, which decomposes the surface into 24 conformally 
equivalent triangles, of angles ^, ^, ^. Therefore an algebraic equation for the 
Lawson surface is = a;^ — 1, with six branch points at the roots of unity. The 
correspondance between the points in the square picture of the surface and the 
double sheeted cover of the complex plane is done in Fig. [3l In particular the 
center of the six squares are sent to the branch points, the vertices are sent to 
the two copies of (black and dark gray) and oo (white and light gray), the 
square are sent to double sheeted two gons corresponding to a sextant. 

Similarly to Sec. l4.2l we compute the period matrices Jlex and Ojn for different 
resolutions utilizing weights imposed by the extrinsic and intrinsic metric and 
compare the results with our conjectured period matrix for the Lawson surface 



^vertices 






1162 


1.68- 10-^ 


1.68-10-'' 


2498 


3.01 • 10-^ 


3.20- 10-^ 


10090 


8.55-10-^ 


8.56 • 10-'^ 



Our first observation is that the matrices ilex and flm almost coincide. Hence 
the method for computing the p seems to have only little influence on this result 
(compare also Sec. 14. 2p . Further we sec that figures of the higher resolution 
surface, i.e. with 2498 and 10090 vertices ar worse than the coarsest one with 
1162 vertices. The mesh beautifying process was most successful on the coarsest 
triangulation of the Lawson surface (Fig. [2|). The quality of the mesh has 
a significant impact on the accuracy of our computation: One can sec that 
the triangles on the coarsest example are of even shapes with comparable side 
lengths, while the finer resolution contains thin triangles with small angles. The 
convergence speed proven in [19] is governed by this smallest angles, accounting 
for the poor result. Therefore for this method to be applicable, the data should 
be well suited and it is not enough to have very refined data if the triangles 
themselves are not of a good shape. A further impediment to the method is the 
fact that the triangulation should be Delaunay. If it is not, it can be repaired 
by the algorithm described in [15j . 
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Figure 3: The Lawson surface is conformally equivalent to a surface made of 
squares. 
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